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ABSTRACT \ 

Quantized vortices can occur around nodal points in wavefunctions . 
This fact, discovered by Dirac (1931) but little noted since, is 
rederived here and examples are discussed. The derivation depends only 
on the wavefunction being single valued, continuous, and having 
continuous first derivatives. Since the derivation does not depend 
upon the dynamical equations, the quantized vortices are expected to 
occur for. many types of waves (i.e. electromagnetic, acoustic, etc.). 
Such vortices have appeared in the calculations (McCullough and Wyatt, 
Kuppermann) of the H + H 2 molecular collisions and play a role in 
the chemical kinetics. In a companion paper, it is shown that 
quantized vortices, occur when optical waves are internally reflected 
from the face of a prism or particle beams are reflected from potential 
energy barriers. 
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whirlpool effect" when they integrated the quantum mechanical equations 
of motion for the col linear collision of a hydrogen atom and a hydrogen 
molecule. This was a vortex in the flux around a node in the 
wavefunction of the composite system moving on a potential energy 
surface. From a classical point of view, it is surprising that the 
flux might go around in closed loops when the potential surface does 
not have a dimple. Since the vortex occurred near the saddle point in' 
the potential surface, the whirlpool effect may. be important in the 
chemical kinetics. (The vortex produces multiple crossings of the 
activated state and thereby lessens the steric factor in the theory of 
absolute reaction rates.) Recent unpublished calculations of 
Kuppermann, and also Wyatt, are said to show interesting variations 
in the behaviour of the vortex with energy, especially near resonances. 

In the preceding paper, 2 vortices are found surrounding nodes in 
the wavefunction describing an optical beam internally reflected at a 
prism face. Or, with another interpretation of the mathematical 
treatment, vortices are found surrounding nodes in the wavefunction 
describing the reflection of particles impinging upon a square 
potential barrier. Clearly, the occurrence of yprtices is a very 
general phenomenon. These quantum whirlpools would not be evident 
from classical trajectory calculations since they result from the 
interference of waves. 

Here we present a rederivation of a fact discovered in 1931 by 
Dirac: 3 Vortices are to be expected around the nodal points of 
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wavsfunctions and t.^ese vortices are quant i zed. We use the terrr! 
"wave-function 11 to denote the spin-free component of a wavefunction. 
The proof is elementary and depends only upon the requirement that 
the wavefunction be single-valued and continuous and its first 
derivatives be continuous. There is no requirement that the 
wavefunction satisfy any particular dynamical equations. Thus, the 
same type of quantized vortices should occur in electromagnetic, 
acoustic, hydrodynamic, and other types of wave mechanics. Our 
formulation is expressed in terms of Schrodinger quantum mechanics, 
but the equations and examples apply equally well to these other 
types of waves provided that the meanings and values of the constants 
f) and m are changed appropriately. 

Our discussion of nodes and vortices also applies to superfluids 
or superconductors which are characterized by a "macroscopic 
wavefunction" or amplitude Y(r,t) . The density p and the local 
mean velocity v are defined in terms of ¥(r,t) by the same 
equations as in Schrodinger wave mechanics. However, the superfluids 
or superconductor ¥ differs from the Schrodinger wavefunction in 
not satisfying a linear equation (i.e. a linear combination of two 
solutions and v 2 is not itself a solution). For superfluids 

or superconductors, the Madelung 4 type of hydrodynamical formulation 
of the quantum mechanical equations is convenient, whereas, for 
Schrodinger wave mechanics, the hydrodynamical formulation has had 
only limited usefulness up to the present time. 
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I. DERIVATION 

As in the preceding paper, 2 we express an arbitrary Schrodinger 
wavefunction in the form y(r,t) - p^ 2 exp ( i S/ft ) where both p^ 2 
and S are real. Then the probability density is p = y*y and the 
phase of y is S/fT . The number flux is 

J = (Hf/2m)[yvy* - y*vy] , (1) 

where m is the mass of the particles. The Madelung-Landau-London 
definition of the local mean velocity is 

v - J/p = (1/m) im[vy/y] = vS/m . (2) 

The fluid flow is irrotational around points where y is not equal 
to zero, since there v * v = v x vS/m - 0 . However, around 
nodes in the v/avef unction, v becomes indeterminate since both 
J - 0 and p = 0 . It is this indeterminancy which leads to the 
possibility of vortex motion around nodes. 

Both y and vy are required to be continuous and single 
valued everywhere except that at points where the potential energy 
is infinite, y remains continuous but vy may be discontinuous. 

It follows that we can, and do, choose so that it is everywhere 

y 9 

continuous and vp z so that it is continuous except at points where 
the potential energy is infinite. (Mote that this choice implies a 
restriction on the value of S since y remains unchanged if the 
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sign of p z is changed at the same time that ?r is added to the 
value of S/H .) When y = 0 , neither S/ff nor vS are well 
defined. At all points where ¥ does not vanish, vS is well 
defined; and at all points where 'a is not zero, S/ff is defined 
up to an integer multiple of 2n . 

The existence of vortices can be explained in the following 
manner: If Sp is the principal value of S and n is an integer 

then Sp/fi and (Sp/fi) + 2-rrn are indis tinguishable since 
exp[i$p/H] = exp [i (Sp/fi) + 27rni] . Because of this lack of uniqueness 
of S/H , the line integral of v between two points r x and r 2 is 

r r 2 r*2 

v*dr = m 1 vS-dr * m" 1 [S p (r 2 ) - S^r^] + 2Tffi/m)n . (2) 

The value of n depends upon the path which is followed in going 
from ri to r 2 . In following along any particular path (not 
passing through a node) vs is well defined so that the value of the 
line integral is unique. If two paths between ri and £2 lie 
close together and neither passes through a node, then the value of 
the corresponding line integrals must be the same. One might expect 
that these two line integrals might differ slightly since vs can be 
expected to have slightly different values at comparable points along 
the two paths. However, according to Eq. (3), the values of these 
two line integrals can differ only by integer multiples of (2nfi/m) 
and, hence, cannot differ slightly. Thus, we reach the conclusion 
that any two line integrals between the same end points have the same 
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value provided that one path can be deformed into the other without 
intersecting any region where ¥ = 0 . Furthermore, if the end points 
are coincident, so that the paths are closed loops, the location of 
the end point is irrelevant and Eq. (3) becomes 

0 V'dr = (2Tffi/m)n . (4) 

For any two closed loops, the value of the line integral is the same 
provided that the loops can be deformed into one another without 
intersecting any point where ¥ = 0 . For example, if a closed loop 
does not enclose any nodal points, then the closed loop can be 
deformed by shrinking it into a circle of radius zero without changing 
the value of the line integral. In such a case, the value of the line 
integral must be zero, or n = 0 , in agreement with our previous 
conclusion that the fluid flow is irrotational except in the vicinity 
of a nodal point. Conversely, there can be no angular velocity around 
a point for which v does not vanish. Atomic wavef unctions provide a 
good example: Only S-states have p nonvanishing at the nucleus. 

If the closed loop encloses a node in the wavefunction at the 
point r 0 » then the value of |n| depends upon the behaviour of the 
wavefunction in the vicinity of the node. Let us take the line 
integral over a circular path of radius q with r 0 as center. The 
circular path is taken in an xy-plane where x = q cos<f> and 
y = q sin<J> . The values of ¥ and S at a point on the path 
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corresponding to the angle 4> can be designated as and 5(<j>) 

respecti vely . Here 

S ( 4>)/1T = n ' 7T + tanp" 1 {Im[¥(^)]/Re['F(^)]} s (5) 

where n' is an integer and tanp" 1 is the principle value of the 
arc tangent. The tan[S($)/tf] becomes discontinuous at a set of 
points <j> a <f>k where 

Re[ l p(<j> k )]/ImC'i f {0k) 3 = 0 . (6) 

Let Sfc be the sign of tan[S ^ - e)^T] where e is a small angle, 
then if S($) is to remain continuous, each time that <j> passes 
through a ^ for which is positive, n ! increases by one; 
and every time it passes through a for which s^ is negative, 
n' decreases by one. Thus, from Eqs. (3) and (4), 

ffl v-dr = (-n/m) [S (2-rr) - S(*)] = (irfi/m) l s k 1 = ( 2flTf/m) n . (7) 

Here the summation is over all of the 4 ^ between 4 = 0 and 2-n . 
Since v is single-valued, there is an even number of <j»^'s so 
that n is an integer. 

The existence of quantized vortices around nodes is well known 
for analytical functions in the complex plane. However, it is 
surprising that the same sort of behaviour also applies to complex 
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wavef unctions in an arbitrary number of dimensions D provided that 
they are single valued, continuous, and have continuous first 
derivatives. As we have seen, there is associated with each node of 
'F a magnitude of circulation |n| corresponding to the circulation 
around a simple loop which encloses this particular node and no other 
node. The sign of the circulation depends upon the direction in which 
one progresses around the loop. Of course, the circulation |n| of 
this node can be nonvanishing only if the node is topologically 
equivalent to a (D-2)-dimensional subspace without a boundary. 
Otherwise, the wavefunction could not satisfy its single-value and 
continuity requirements. For instance, in three dimensions, a node 
with circulation must be a ring, an endless ("vortex") line, or else 
a ribbon or rod which is topologically equivalent (can be shrunk) to 
an endless line--a line with end points or a single point would not 
suffice. 

The wavefunction can be expanded in a many-dimensional Taylor's 
series about a nodal point ro * Then, if the node is of the g-th 
order so that the first nonvanishing term involves the g-th order 
of the coordinates, it is easy to prove (by examining the possible 
number of the ^'s ) that |n| must be less than or equal to g . 

The first order nodes occur most frequently. Near a first order 
nodal point r 0 , the wavefunction can be expanded in the form 




[a+ig]-(r-r 0 ) + [y+is]:(r-r 0 )(r-r 0 ) + ... s 


( 8 ) 
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where the a, 3 , y, ^nd 6 are real and constant. In order that the 

node be of the first order, either a or g must be nonvanishing. 

The second order term (involving y and <5 ) and the higher order 

terms then do not affect the value of n . If a and 6 are 

parallel so that 6 = ca , this would be an example of a 

(D-l ) -dimensional node since, to the first approximation, the only 

nonvanishing configurations of r - r 0 lie along a line in the 

direction of a . Since tan(S/H) = c , it is obvious that n = 0 

which is confirmation of the theorem that a (D-l )-dirnensional node 

is topologically unsuited to support a vortex. The same argument 

would apply to the case where either a or 3 were zero. If a, 

and 3 are not parallel, we can choose the x-axis to be in the 

direction of a and choose the y-axis so that 6 lies in the 

xy-plane. Such a node is (D-2)-dimensional since the only 

nonvanishing configurations of (r - r 0 ) lie in the plane defined 

by a and 3 . For this case, tan(S/tT) = (3 ¥ /ct ) + (g / a }tan<j> . 

xx y x 

Thus, there are two values of for which tan(SAti) becomes 
discontinuous, <Pi - v/2. and <f >2 = 3tt/ 2 , and S} = S 2 so that 
either n = 1 or n = -1 . If vy is parallel to V¥* in the 
vicinity of a node, then that node cannot support a circulation. A 
case of this nature is ¥ = y exp(ikx) . Here S/H « kq cost so 
that tan(S/ft) has no discontinuities and n = 0 . 

It is easy to show that in the presence of a magnetic field with 
a vector potential A , the extra term, -(e/mc)A , in v does not 
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contribute to the value of the line integral provided that the radius 
q is required to approach zero* 

1 1 • EXAMPLE: AM INFINITE FIELD OF REGULARLY SPACED VORTICES 

An infinite number of vortices with m - 1 or m ~ --1 appear 
at regular spacings when two plane waves, H'j and y 2 * having the 
same energy E = .k 2 fi 2 /2y but corresponding to motions in different 
directions, are reflected from a plane wall so that their reflected 
wave-functions are '?{ and y 2 . Thus, the total wavef unction for 
the system is ¥ = Y 1 +y 2 + 4 , { + 4 f 2 . For simplicity, let the plane 
wall be the surface x - 0 and the incoming wavefunctions be 


4*! “ i exp[-ik(x cos^ + y sin^)] 


and 


( 9 ) 


t 2 " i exp[-ik(x cos <j> 2 + y sin<f> 2 )3 » 

where j and <£ 2 are constants. The reflected wavefunctions are 
then 


and 


= -i exp[ik(x cos^ - y sin^)] 
^2 " exp[ik(x ces<j> 2 - y s 1 n 2 ) J - 


( 10 ) 


For notational convenience, let 5 = k cos^ , b = k $in<f> , 
a = k cos<p z } and 3 = k sin<£ 2 . Then, 
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Y = 2 sin (ax) exp (-i by) +2sin(ax) exp(-igy) (11) 

so that there are four conditions under which y = 0 , 

(1) . cos[(b-g)y 0 ] = 1 and si n[(a+a)x 0 /2] = 0 so that 

x 0 - 2 tt j / (a+a) , y 0 = 2Trk/(b-g) ; 

(2) . cos[(b-g)y 0 ] = -1 and sin[(a-a)x 0 /2] = 0 so that 

x 0 s 2 ttj/ ( a-ot) , y 0 = tt ( 1 -f-2k)/ ( b-3 ) ; 

(3) . cos[(b-g)y 0 ] = 1 and cos[(a-a)x 0 /2] = 0 so that 

*o = 7r(l+2j)/(a-a) , y 0 = 27rk/(b-g) ; 

(4) . cos[(b-$)y 0 ] = -1 and cos [(a+a)x 0 /2] = 0 so that 

x 0 a "(l+2j)/(a+a) , y 0 - 7T ( 1 +2k )/(b-6 ) , 

where the constants j and k are either zero or a positive cr 
negative integer. 
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